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Abstract 

We have examined the electrical transport in disordered molecular solids. It has been found 
that mobility is a function of electric field and temperature due to hopping conduction. Several 
theoretical models for charge transport in disordered solids have been debated over the role of 
spatial and energetic correlation in these systems and such correlations have been recently shown 
to explain the universal electric field dependence of mobility. We have compared and evaluated 
the applicability of different theoretically proposed models using very simple experimental results 
and based on our extensive analysis, we have found that correlation is important to explain the 
electrical transport in these systems. 
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Recently, considerable debate 



arrier transport in 

disordered organic molecular solids(DOMS) because of two reasons: first, their application 
in display devices and second, the fundamental understanding of charge transport in these 
materials. There are several distinguishing features of these solids, (i) they are composed of 
organic molecules, held together loosely by weak van der Waals type intermolecular coupling 
while the intramolecular coupling is strong; (ii) the absence of long range order in these 
disordered materials lead to the localization of the electronic wave function and the formation 
of a broad Gaussian density of states(GDOS) and the most important one is (iii) the carrier 
mobility /x exhibits a nearly universal Poole-Frenkel(PF) behaviorjj, ^ 



/i(£;)=/i(0)exp(7v^) (1) 

where /i(0) and 7 are temperature dependent quantities, known as the zero field mobility and 
the field activation of the mobility, respectively. In DOMS, hopping among the molecular 
sites having comparable energies, describes the transport of charge carriers through the 
GDOS of highest occupied molecular orbital(HOMO) and/or lowest unoccupied molecular 
orbital(LUMO). 

There are several models to explain the universal feature of carrier mobility /i(Eq.l) ob- 
served in wide class of organic solids. Gillj^ first attempted to describe the experimental 
data with the following empirical form of zero field mobility yu(0) and field activation of mo- 
bility 7, /i(0) = /ioexp[— A/ {ksT)] and 7 = B[l/{kBT) — 1/(/cbTo)], where /iq is the zero field 
temperature independent mobility, A is the zero-field activation energy, and B and Tq are 
parameters of the model. There are several conceptual problems with this phenomenologi- 
cal description and the most important are (i) the empirical Gill's formula lacks theoretical 
justification and the significance of different parameters(A, B and To) are not well under- 
stood in the context of the physical properties of the system and (ii) the empirical formula 
for 7 predicts the negative field dependence for T > Tq, which deviates from the universal 
feature(PF) of carrier transport in these materials. It is now recognized that the spatial 
fluctuations in the potential energy of the charge carriers in DOMS results the PF behav- 
io. ot ca.„e. o^obmtyQQQ. There a.e two co,npet.ng o^odel, debating ovejthe to.e of 
correlation on the spatial fluctuations in these systems. Bassler and co-workers [2I| proposed 
the uncorrelated Gaussian disorder model(UGDM) and provided support to PF behavior 
of mobility using Monte Carlo simulation. The UGDM describes the carrier transport as 
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a biased random walk among the dopant molecules with Gaussian-distributed random site 
energies. UGDM lead to following form of /i(0) and 7 

/i(0) = /ioexp |- ,7 = C{a^ - S) (2) 

where a = cr/kT, a is the width of the Gaussian distribution of the hopping sites energies, 
and C and S are parameters of the model. In this case, /i(0) and 7 have depen- 
dence, instead of 1/T dependence as in case of Gill's empirical formula. Although UGDM 

Q 

explains |7| some features of experimental data and provides support for PF behavior of 
carrier mobility, several discrepancies emerge with uncorrelated description of Gaussian dis- 
order model, which will be discussed later. The most important criticism against UGDM 
is its' inability to reproduce the PF behavior over wider range of electric field. Garstein and 
Conwellj^ first showed that a spatially correlated potential is required for the description 



of PF behavior of mobility for wider range of electric field. Dunlop and co-workers [ 
have shown analytically in ID and numerically in 3D that the interaction of charge carriers 
with permanent dipoles located on either dopant or host molecules give rise to PF behav- 
ior of mobility. Essentially, this correlated Gaussian disorder model(CGDM) is based on 
long-range correlation between charge carriers and the molecular electric dipole. They have 
shown that this interaction gives rise to random potential energy landscape with long-range 
spatial correlations < f/(0)f/(r) >~ a'^a/r, where a is the minimal charge-dipole separation 
or the lattice constant of the molecular solid. CGDM gives the following form of /i(0) and 7 

m = ^oexp |- (^^)'} ,7 = A {a^^ - v) (^)'^' (3) 

where A and F are parameters of the model. Here, /i(0) has a dependence (same as 

in UGDM) but 7 has a dependence. In this letter, we present the measurements 

on charge carrier mobility fi and its field activation 7 as a function of temperature and 
electric field. The experimental results are analyzed within UGDM and CGDM. It has 
been found that correlated transport model CGDM successfully explains the experimental 
data emphasizing the important role of correlation in explaining the transport mechanism 
in disordered organic solids. 

We have chosen two different metal-phthalocyanines(MePc), Copper- 
phthalocyanine(CuPc) and Zinc-phthalocyanine(ZnPc) for our investigations because 
of their high chemical stability, ease to prepare the devices and reproducibility of ex- 
perimental data, which is a major problem with most other organic solids due to their 



degradation with time. Details about the MePc-based single layer devices are given in 
Ref. 11 1 and Ref.^^. In this letter, we report the experimental investigation on charge car- 
rier transport in hole only devices based on metal/MePc/metal structures. In this case, by 
properly choosing contacting metals, current injection and transport due to holes have been 
investigated. Indium tin oxide(ITO) has work function(4.75eV) very close to the ionization 
potential(4.8eV) of MePcQ, resulting the ITO/MePc interface as an Ohmic contact. 
Fig.l shows the current- voltage( J- V) characteristics of ITO/CuPc/Al and ITO/CuPc/Cu 
at room temperature. The experiment consisted of two steps. First: the current due to 
hole injection from ITO was measured by supplying it to positive bias and second: the 
current due to hole injection from Cu and Al was measured by reversing the polarity of 
the bias voltage(i.e., biasing Al and Cu electrodes positively). It is clear from Fig.l that 
J-V characteristics in case of ITO/MePc/Cu display almost symmetric behavior in both 
cases(hole injection either from ITO, or from Cu electrodes). This is because in both cases 
there is small energy barriers of O.OSeV in case of ITO/MePc and 0.2eV in case of Cu/MePc 
interface, giving rise to space charge limited(SCL) bulk current when either ITO or Cu is 
positively biased. It has been shown , for Schottky energy barrier(SEB) less than 

MUM. But, in case of ITO/MePc/Al 



about 0.3-0.4eV, the current flow is due to SCL 
devices, J-V characteristics displays asymmetric and rectification-like behavior. Current 
density increases by almost five order of magnitude by reversing the polarity of the bias. 
As discussed before, when ITO is positively biased, current is due to SCL and when Al 
is positively biased, current is due to injection limited and reduced by several order of 
magnitude due to the existence of SEE of 0.6eV at MePc/Al contact. It has been shown by 



us 



llj that this injection limited current J^j from a metal electrode to disordered organic 
semiconductor can be described by the relation, Jinj oc -ip"^ nexp{—(l)B / kBT)exp{f^^'^) , ip is 
a slow varying function of E, (ps is SEB and / = e^E/2TTe{kBT)'^, proposed by Shen et. 
ai.H. 

Inset of Fig.l shows the temperature dependent SCL current in ITO/MePc/Al structure 
made out of 200nm thick active organic layer of CuPc. At low bias, J-V characteristic 
follow the slope of 2 in log(J)-log(V) plot, but as the bias increases the slope increases 

r I 

gradually to higher value. A numerical workout [13] has been used by solving Poission's 
equation, dE{x)/dx = ep{x)/e and J{x) = ep{x)fi{x)E{x) simultaneously, and assuming 
PF form(Eq.l) of the carrier mobility. Here, e is the dielectric constant of MePc and p{x) 
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is the hole density at position x and the boundary condition is taken at the Ohmic con- 
tact (ITO/MePc interface) with hole density(at x = 0) of 2.5 x lO^^cm"^, which is the 
density of states of HOMO in MePc. The simulated results are shown as solid lines in in- 
set of Fig.l at different temperatures. The excellent agreement of the simulated and the 
experimental data confirms PF behavior of carrier mobility and the thermally activated 
hopping transport. Fig. 2 shows the temperature dependence of conductivity <j(T) for CuPc. 
Linear dependence of log{a) on validates the variable range hoppingfVRH) [isj] . 

cr(T) = croexp{ — (Ti/T)^/^}, where cxo is the temperature independent conductivity and 
Ti is the characteristic temperature. The inset shows the dependence of ctq at different elec- 
tric field and linear dependence of log{aQ) on ^/E signifies the PF behavior in charge carrier 
transport conclusively. Hence the experimental data presented in Fig.l and Fig.2 establish 
the carrier transport process, which is thermally activated hopping with field dependent 
mobility. 

Now, we compare our experimental data with the uncorrelated(UGDM) and corre- 
lated(CGDM) theoretical models proposing the stretched exponential electric field depen- 
dence of /i. At every temperature, /i(0) and 7 are obtained from the experimental data and 
plotted in Fig. 3 and Fig.4. Fig. 3 shows the temperature dependence of the /i(0) and 7 ac- 
cording to UGDM and the apparent linear dependence of ln/i(0) on and 7 on can 
be described by Eq.(2). Values of temperature independent mobility /xq, width of the energy 
spread a, and the constant C and S are determined from the ln/i(0) vs. and 7 vs. 

plots, and given in Table I. Main feature of UGDM is the non-Arrehenious behavior of 
the temperature dependence of mobility, which is the consequence of hopping conduction in 
GDOS. But, there are several discrepancies in this description, (i) fitting is not as good as in 
case of CGDM, which will be discussed next; (ii) the intercept of 7 vs. l/T^(shown in Fig. 3) 
gives positive value resulting negative value of parameter E(given in Table I), but according 
to UGDM(Eq.2), it should be negative, which is a fundamental problem with UGDM, and 
(iii) it has been observed that there is no specific trend in the values of C and S, which is 
a scaling parameter in the model and cannot be linked to the physical parameters of the 
system. 

Fig.4 shows the temperature dependence of /i(0) and 7 according to CGDM, described by 
Eq.(3). Values of the parameters fio, a, a and P are determined from the straight line plots 
of ln/i(0) vs. and 7 vs. and given in Table I. Following observations have been 
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made regarding the CGDM description of the experimental data, (i) hnear dependence of 
ln/i(0) on and 7 on is evident in Fig.4 and straight hne fit is excellent and better 

than that in UGDM; (ii) for all samples the intercept of 7 vs. gives negative values, 

according to CGDM(Eq.3); (iii) the experimental value of F are 1.2 in CuPc and 1.4 in ZnPc, 
which are very close to the predicted value 1.97|2J; (iv) the width of the Gaussian distribution 
cr is lOOmeV in CuPc and 120meV in Zn-Pc and similar values for a are commonly observed 
in other molecular solids; (v) the lattice constant (a) has been found to be 19 A in CuPc 

n 

and 17^4 in ZnPc, are in excellent agreement with the reported[19|| experimental values 
determined by X-Ray diffraction and (vi) MePc molecules are symmetric, so the dipole 
moment p should be zero or very small, but the symmetry of the individual molecule may 
be lowered in solid state, so, in any case p should be very small, which is corroborated by the 
estimated small value of p(=0.3D) for CuPc using the relation a = 2.35ep/ea^ , according to 

cgdmQ. 

The UCDM and CGDM share the common feature regarding the temperature dependence 
of yu(0) due to similar distribution of hopping sites energies, but the temperature dependence 
of field activation parameter 7 decides critically the importance of spatial correlation and 
decides the mechanism of carrier transport in organic molecular solids. Though CGDM 
explains the experimental data successfully, there are some limitations with this model. The 
positional disorder is not included in the present form of CGDM, but large values of disorder 
parameters a and F in NPPDA ^, A1q[20| and MePc signify the importance of the positional 
disorder. In CGDM, it has been shown that long-range interaction between charge carriers 
and permanent dipole moments of doped molecules in polymers and host molecules leads to 
spatial correlation. However, it has been pointed out by Yu et. al. 2l| that the mechanism 
responsible for PF behavior in different conjugated polymers and molecules cannot be due to 
only charge-dipole interaction, because PF behavior has been universally observed in several 
doped and undoped polymers and molecules with or without permanent dipole moment. 
Hence, in addition to charge-dipole interaction there may be another mechanism responsible 
for spatial correlation, which is a fundamental requirement for PF behavior. Yu et. ai.js^ 
have shown using first principle quantum chemical calculation that the thermal fluctuations 
in the molecular geometry can lead to spatial correlation. It has been shown that the 
primary restoring force for the thermal fluctuation is steric and intermolecular, which lead to 
spatially correlated fluctuation in the energies of the localized states. Further experimental 
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investigations are in progress to identify the exact origin of spatial correlation in molecules 
with and without dipole moments. 

In summary, we have studied the mechanism of carrier transport in disordered molecular 
solids by means of simple experiments. We have (i) shown that the carriers conducts through 
hopping; (ii) given evidences for the universal electric field dependence(/n(/x) cx VE) of car- 
rier mobility and (iii) presented a comparison of the transport properties with uncorrelated 
and correlated disorder model. We have shown that E and T dependence of experimental 
data can be fitted and described within the correlated Gaussian disorder model, signifying 
the important role of both energetic and spatial disorder on the carrier transport in these 
materials. The origin of correlation has been discussed. Wc hope present analysis would 
lead to further experimental and theoretical investigations focusing the exact origin and role 
of correlation on transport properties of disordered molecular solids. 

Table.I. Temperature independent zero field mobility /i(0), width of energy spread a(ai in 
UGDM and an in CGDM) and different scaling parameters C, E in UGDM and a, F in 
CGDM for CuPc and ZnPc are given. 
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Figure caption. 

Fig.l. Room temperature J-V characteristics for single layer hole only CuPc device with 
thickness of 200nm. Empty circles represent the data when ITO is positively biased 
and empty squares represent the data when ITO is negatively biased. Inset shows the 
J-V characteristics of bulk limited current in ITO/CuPc/Al at different temperature 
staring at 320K and then at the interval of 40K. Empty circles show experimental data 
and solid lines represent the simulated data. 

Fig.2. Temperature dependence of conductivity of CuPc. Empty circles are experimental 
data and solid hne is fit with straight line. Inset shows the electric field dependence 
of temperature independent conductivity. Solid line is fit with PF dependence of 
mobility. 

Fig. 3. Temperature dependence of the field activation parameter 7 and zero field mobil- 
ity /i(0) (inset) for CuPc obtained from the theoretical fit to the experimental data, 
according to the UGDM(Eq.2). 

Fig.4. Temperature dependence of the field activation parameter 7 and zero field mobil- 
ity /x(0) (inset) for CuPc obtained from the theoretical fit to the experimental data, 
according to the CGDM(Eq.3). 
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